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We carry out a high-precision Monte Carlo simulation of the two-dimensional 0(3)-invariant fr-model at cor- 
relation lengths { up to ~ 10 . Our work employs a new and powerful method for extrapolating finite- volume 
Monte Carlo data to infinite volume, based on finite-size-scaling theory. We compare the extrapolated data to the 
renormalization-group predictions. The deviation from asymptotic scaling, which is Si 25% at { ~ 10 , decreases 
to w4% at { ~ 10^ 



We study the lattice cr-model taking values 
in the unit sphere 5'^""'^ C R , with nearest- 
neighbor action 'H{(t) = —(3 '}2'^x ■ o'y Perturba- 
tive renormalization-group computations through 
three loops [1,2] predict that the exponential cor- 
relation length (inverse mass gap) ^(^^p)^ the 
second-moment correlation length S} \ and the 
susceptibility x behave for # = 3 (in infinite vol- 
ume) as 
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as /3 ^ cxD. The nonperturbative constant Cf(e^p) 
has been computed recently using the thermody- 
namic Bethe Ansatz [3]: 
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The remaining nonperturbative constants are 
known analytically only at large N [4] : 
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A high-precision Monte Carlo study [5] yields 
C^(2)/C^(e.p) = 0.9993 ± 0.0006 for # = 3, in 
good agreement with (4). Previous studies up to 
^ ~ 100 agree with these predictions to within 
about 20-25% [6]. 

In order to extrapolate finite-volume Monte 
Carlo data to infinite volume, we used a novel 
method [7] based on the finite-size-scaling Ansatz 
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which is correct up to terms of order ^~^ and 
L""^; here O is any long-distance observable, s 
is a fixed scale factor (usually s = 2), _L is the 
linear lattice size, Fq is a universal function, and 
w is a correction-to-scaling exponent. For similar 
extrapolation methods, see [8,9]. 

Details of our simulation and of the extrapola- 
tion process can be found in Sokal's talk at this 
conference and in [10]. Our preferred fit is shown 
in Figure 1, where we compare also with the per- 
turbative prediction 
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valid for x '^ I. 
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Figure 1. ^{p,2L)/^{p,L) versus ^{(],L)/L. Dif- 
ferent symbols indicate different lattice sizes L. 
Dashed curve is the perturbative prediction (7). 



The extrapolated values ^Jx) from different lat- 
tice sizes at the same /3 are consistent within sta- 
tistical errors: only one of the 24 /3 values has a 
X^ too large at the 5% level; and summing all /3 
values we have x^ = 86.56 (106 DF, level = 92%). 

In Figure 2 (points -|- and x) we plot 

^''^lestimate divided by the two-loop and three- 
loop predictions (l)-(4). The discrepancy from 
three-loop asymptotic scaling, which is Ri 16% at 
l3 = 2.0 (C ^ 200), decreases to ^ 4% at /3 = 3.0 
{£, K, 10®). This is roughly consistent with the ex- 
pected l//3^ corrections. Notice that the points 
in the curve fluctuate less than what one would 
expect on the basis of the reported error bars; this 
is due to a strong statistical correlation of the es- 
timates for the points with higher values of (3. 
Probably also the slight bump at 2.3 ^ /3 ^ 2.6 
is spurious, arising from correlated statistical or 
systematic errors. 

We can also try an "improved expansion pa- 
rameter" [11,2,13] based on the energy E = (<to • 
<Ti). First we invert the perturbative expansion 
[12,2] 
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Figure 2. S,Z\estimatel^%heor ^ersus /3. Error 
bars are one standard deviation (statistical error 
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only). There are four versions of ^^ 
dard perturbation theory in 1//3 gives points 
(2-loop) and x (3-loop); "improved" perturba- 
tion theory in 1 — £" gives points □ (2-loop) and 
O (3-loop). 



and substitute into (1); this gives a prediction 
for ^ as a function of 1 — £". For E we use the 
value measured on the largest lattice; the statisti- 
cal errors and finite-size corrections on E are less 
than 5 x 10~®, and therefore induce a negligible 
error (less than 0.5%) on the predicted ^. The 
corresponding observed/predicted ratios are also 
shown in Figure 2 (points □ and O). The "im- 
proved" 3-loop prediction is in excellent agree- 
ment with the data. 

In Figure 3 we report the ratio 
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where by the suffix theor we denote the one-, two- 
or three-loop prediction either for the standard 
perturbation theory in 1/(3 or for the "improved" 
one in l — E. The curves appear to become fiat for 
increasing /3 and to converge to the same value Ri 
10.8, thus providing an estimate for the universal 
ratio C;^/C?(2). 
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Figure 3. R^ versus /3. There are six versions 
of perturbation theory: the standard one in 1//3 
gives points + (1-loop), x (2-loop) and t±t (3- 
loop); "improved" perturbation theory in 1 — £" 
gives points □ (1-loop), O (2-loop) and tf (3- 
loop). Error bars are one standard deviation and 
are reported only for one series of prediction be- 
cause they are always the same at fixed /3. 



Let us summarize the conceptual basis of our 
analysis. The main assumption is that if the 
Ansatz (6) with a given function F^ is well sat- 
isfied by our data for Lmin ^ L < 256 and 
1.65 < /3 < 3, then it will continue to be well 
satisfied for L > 256 and for /3 > 3. Obviously 
this assumption could fail, e.g. if [14] at some 
large correlation length (> 10^) the model crosses 
over to a new universality class associated with a 
finite- /3 critical point. In this respect our work is 
subject to the same caveats as any other Monte 
Carlo work on a finite lattice. However, it should 
be emphasized that our approach does not as- 
sume asymptotic scaling [eq. (1)], as /3 plays no 
role in our extrapolation method. Thus, we can 
make an unbiased test of asymptotic scaling. The 
fact that we confirm (1) with the correct nonper- 
turhattve constant (3)/(4) is, we believe, good ev- 
idence in favor of the asymptotic-freedom picture. 
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